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ABSTRACT: In this paper, we introduce an approach to increase the integration rate of elements of a 
conventional static comparator. Framework, the approach we consider a heterostructure with special 
configuration. Several specific areas of the heterostructure should be doped by diffusion or ion 


implantation. Annealing of dopant and/or radiation defects should be optimized. 
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1. INTRODUCTION 

An actual and intensively solving problem of 
solid state electronics is increasing of integration 
rate of elements of integrated circuits (p-n- 
junctions) [1-8]. Increasing the integration rate 
leads to necessity to decrease their dimensions. To 
decrease the dimensions are using several 
approaches. These are widely using in the laser 
and microwave types of annealing of infused 
dopants. These types of annealing are also widely 
used for annealing of radiation defects, generated 
during ion implantation [9-17]. Using the 
approaches gives a_ possibility to increase 
integration rate of elements of integrated circuits 
through inhomogeneity of technological 
parameters due to generating inhomogenous 
distribution of temperature. In this situation one 
can obtain decreasing dimensions of elements of 
integrated circuits [18] with account Arrhenius 
law [1,3]. Another approach to manufacture 
elements of integrated circuits with smaller 
dimensions is doping of heterostructure by 
diffusion or ion implantation [1-3]. However in 
this case optimization of dopant and/or radiation 
defects is required [18]. 


In this paper, we consider a heterostructure, this 
consist of a substrate and several epitaxial layers. 
Some sections have been manufactured in the 
epitaxial layers. Further we consider doping of 
these sections by using diffusion or ion 
implantation. The doping gives a possibility to 
manufacture field-effect transistors framework a 
cascaded-inverter circuit (See Figs. 1). The 
manufacturing gives a possibility to increase 
density of elements of the conventional static 
comparator [4]. After considered doping dopant 
and/or radiation defects should be annealed. 


Framework, in the paper we have analyzed 
dynamics of redistribution of dopant and/or 
radiation defects during their annealing. We 
introduce an approach to decrease dimensions of 
the element. However it is necessary to complicate 
technological process. 


Fig. 1 The considered cascaded-inverter [4] 


METHOD OF SOLUTION 


In this section we determine the _ spatio- 
temporal distributions of concentrations of 
infused and implanted dopants. To determine 
these distributions we calculate appropriate 
solutions of the second Fick's law [1,3,18] 
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Boundary and initial conditions for the equations are 
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The function C(x,y,z,t) describes the spatio- 
temporal distribution of concentration of dopant; 
T is the temperature of annealing; Dc is the dopant 
diffusion coefficient. Value of dopant diffusion 
coefficient could be changed with changing 
materials of heterostructure, with changing 


D.=D,(x,y,z,T)/1+é 


Here the function Dz (x,y,z,T) describes the spatial 
(in heterostructure) and temperature (due to 
Arrhenius law) dependences of _ diffusion 
coefficient of dopant. The function P (x,,z,T) 
describes the limit of solubility of dopant. 
Parameter y €[1,3] describes average quantity of 
charged defects interacted with atom of dopant 
[20]. The function V (x,y,z,t) describes the spatio- 
temporal distribution of concentration of 
radiation vacancies. Parameter V’ describes the 


C’(x,y,z,t) 
P’(x, V2; T) 


x=L, 


temperature of materials (including annealing), 
with changing concentrations of dopant and 
radiation defects. We approximate dependences of 
dopant diffusion coefficient on parameters by the 
following relation with account results in Refs. 
[20-22] 


2 
V ae 3) 


equilibrium distribution of concentration of 
vacancies. The considered concentrational 
dependence of dopant diffusion coefficient has 
been described in details in [20]. It should be 
noted, that using diffusion type of doping did not 
generation radiation defects. In this situation, (= 
¢2= 0. We determine spatio-temporal distributions 
of concentrations of radiation defects by solving 
the following system of equations [21,22] 
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Boundary and initial conditions for these equations are 


dplxy.2t)) _9 dp(xy zt) _,. Geluyzt)) _, Gelnyzt} _» 
Ox Ss Ox Se Oy is Oy = 
Cale. s.est) 0, 0px yz) = 0, p(xyz,0)=fo&%y.2). (5) 


Here o =1,V. The function I (x,y,z,t) describes the 
spatio-temporal distribution of concentration of 
radiation interstitials; D,(x,y,z,T) are the diffusion 
coefficients of point radiation defects; terms 
V2(x,y,z,t) and /2(x,y,z,t) correspond to generation 
divacancies and diinterstitials; kiv(x,y,z,T) is the 
parameter of recombination of point radiation 


defects; kii(x,y,z,T) and kv v(xy,z,T) are the 
parameters of generation of simplest complexes of 
point radiation defects. 


Further, we determine distributions in space and 
time of concentrations of divacancies ®/(x,y,z,t) 
and diinterstitials @(x,y,z,t) by solving the 
following system of equations [21,22] 
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Boundary and initial conditions for these equations are 
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D (x,y,Z,0)=far (xy,Z), Pv (x,y,2,0)=fov (x,y,z). (7) 


Here Do,(x,y,z,T) are the diffusion coefficients of the above complexes of radiation defects; k;(x,y,z,T) and kv 
(x,y,z,T) are the parameters of decay of these complexes. 
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We calculate distributions of concentrations of 
point radiation defects in space and time by 
recently elaborated approach [18]. The approach 
based on transformation of approximations of 
diffusion coefficients in the following form: 
D,(x,y,Z,T)=Do,[1+ & gp{x,y,Z,T)], where Dop are the 
average values of diffusion coefficients, 0<¢,<1, 
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of recombination of point defects and parameters 
of generation of their complexes: 
kiv(X%y,Z,T)=Korv[1+env giv(Xy,Z,T)], kii(xy,zZ,T)=Ko1 
[1+é71 gi(xy,z,T)] and kvv(xy,z,T)=kovy [1+evv 
gvv(x,y,Z,T)], where kopi,.2 are the their average 
values, 0<é,v <1, 0<e1 <1, O<evv<1, | giv(Xy,z,T)|<1, 
| gi(Xy,Z,T)|<1, \gv.v(xy,z,T)|<1. Let us introduce 


lg(x, y,z,T)|<1, pe =1,V. We also used analogous 
transformation of approximations of parameters 


T(x%y,zt)=M~y zt, V(ioyzt)= =Vixoyzt/V', 
Q. =Lk,,, /D,,Dyy a= D,, t/L , x= x/Ls n=y [Ly 6=2/Le. 


the following dimensionless variables: 


o= Deals D,D, , 


The introduction leads to transformation of Eqs.(4) and conditions (5) to the following form 
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We determine solutions of Eqs.(8) with conditions (9) framework recently introduced approach [18], i-e. as 
the power series 


)= £,(1.6,9) 
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Substitution of the series (10) into Eqs.(8) and conditions (9) gives us possibility to obtain equations for 
initial-order approximations of concentration of point defects Tol 1.10.09) and Vino Xo b,9) and 


(x, 1); d, 39), i>1, j 21, k =>1. The equations are presented in 


the Appendix. Solutions of the equations could be obtained by standard Fourier approach [24,25]. The 
solutions are presented in the Appendix. 
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Now we calculate distributions of concentrations of simplest complexes of point radiation defects in 

space and time. To determine the distributions we transform approximations of diffusion coefficients in the 

following form: Do,(x,y,z,T)=Doe,[1+ €&a9 o,(x,,z,T)], where Dow are the average values of diffusion 


coefficients. In this situation the Eqs.(6) could be written as 
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Farther we determine solutions of above equations as the following power series 


® (x, y,z,t)= 2&4, (x y, zt). (11) 


Now we used the series (11) into Eqs.(6) and 
appropriate boundary and initial conditions. The 
using gives the possibility to obtain equations for 
initial-order approximations of concentrations of 
complexes of defects ®,0(x,y,z,t), corrections for 
them @,i(x,y,z,t) (for them i =>1) and boundary and 
initial conditions for them. We remove equations 
and conditions to the Appendix. Solutions of the 
equations have been calculated by standard 
approaches [24,25] and presented in the 
Appendix. 


Now we calculate distribution of concentration of 
dopant in space and time by using the approach, 
which was used for analysis of radiation defects. 
To use the approach we consider following 
transformation of approximation of dopant 
diffusion coefficient: Di(x,y,Z,T)=Do1[1+ 
é1g1(X,y,Z,T)], where Do. is the average value of 
dopant diffusion coefficient, O<ser< 1, 
\g.(x,y,z,T)|<1. Farther we consider solution of 
Eq.(1) as the following series: 


Clu yzt)= Fe ¥e'C, (x y-z.0). 


Using the relation into Eq.(1) and conditions (2) 
leads to obtaining equations for the functions 
Ci(xy,z,t) (i 21, j 21), boundary and _ initial 
conditions for them. The equations are presented 
in the Appendix. Solutions of the equations have 
been calculated by standard approaches (see, for 
example, [24,25]). The solutions are presented in 
the Appendix. 


DISCUSSION 


In this section we analyzed spatio-temporal 
distributions of concentrations of dopants. Figs. 2 
show typical spatial distributions of 
concentrations of dopants in neighborhood of 


We analyzed distributions of concentrations of 
dopant and radiation defects in space and time 
analytically by using the second-order 
approximations on all parameters, which have 
been used in appropriate series. Usually the 
second-order approximations are enough good 
approximations to make qualitative analysis and 
to obtain quantitative results. All analytical results 
have been checked by numerical simulation. 


interfaces of heterostructures. We calculate these 
distributions of concentrations of dopants under 
the following condition: value of dopant diffusion 
coefficient in doped area is larger, than value of 
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dopant diffusion coefficient in nearest areas. In 
this situation one can find increasing of 
compactness of field-effect transistors with 
increasing of homogeneity of distribution of 
concentration of dopant at one time. Changing 
relation between values of dopant diffusion 
coefficients leads to opposite result (see Figs. 3). It 
should be noted, that framework the considered 


L,LyL, 


Leet 
Uae 7 evn2) wna aad vax: (12) 


C(x.) 


Substrate 


Fig. 2a. Dependences of concentration of dopant, 
infused in heterostructure from Figs. 1, on 
coordinate in direction, which is perpendicular to 
interface between epitaxial layer substrate. 
Difference between values of dopant diffusion 
coefficient in layers of heterostructure increases 
with increasing of number of curves. Value of 
dopant diffusion coefficient in the epitaxial layer is 
larger, than value of dopant diffusion coefficient in 
the substrate 
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Fig. 2b. Dependences of concentration of dopant, 
implanted in heterostructure from Figs. 1, on 
coordinate in direction, which is perpendicular to 
interface between epitaxial layer substrate. 
Difference between values of dopant diffusion 
coefficient in layers of heterostructure increases 
with increasing of number of curves. Value of 
dopant diffusion coefficient in the epitaxial layer is 
larger, than value of dopant diffusion coefficient in 
the substrate. Curve 1 corresponds’ to 
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approach one shall optimize annealing of dopant 
and/or radiation defects. To do the optimization 
we used recently introduced criterion [26-34]. The 
optimization based on approximation real 
distribution by step-wise function w (x,y,z) (see 
Figs. 4). Farther the required values of optimal 
annealing time have been calculated by 
minimization the following mean-squared error 


homogenous sample and annealing time © = 
0.0048 (Lx2+L,2+Lz4)/Do. Curve 2 corresponds to 
homogenous sample and annealing time © = 
0.0057 = (Lx*+Ly*+L77)/Do. Curves 3 and 4 
correspond to heterostructure from Figs. 1; 
annealing times © = 0.0048 (Lx2+L,2+L4)/Do and © 
= 0.0057 (Lx4 +Ly?+Lz?) /Do, respectively 


10 4 


C(x.0) 


08 


06 + T T 1 
-a,=-L/3 0 a=2L/3 L 
x 

Fig.3a. Distributions of concentration of dopant, 
infused in average section of epitaxial layer of 
heterostructure from Figs. 1 in direction parallel 
to interface between epitaxial layer and substrate 
of heterostructure. Difference between values of 
dopant diffusion coefficients increases with 
increasing of number of curves. Value of dopant 
diffusion coefficient in this section is smaller, than 
value of dopant diffusion coefficient in nearest 
sections 
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Fig.3b. Calculated distributions of implanted 
dopant in epitaxial layers of heterostructure. Solid 
lines are spatial distributions of implanted dopant 
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in system of two epitaxial layers. Dushed lines are 
spatial distributions of implanted dopant in one 
epitaxial layer. Annealing time increases with 
increasing of number of curves 


C(x,0) 


0 L, 


Fig.4a. Distributions of concentration of infused 
dopant in depth of heterostructure from Fig. 1 for 
different values of annealing time (curves 2-4) and 
idealized step-wise approximation (curve 1). 
Increasing of number of curve corresponds to 
increasing of annealing time 
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Fig.4b. Distributions of concentration of implanted 
dopant in depth of heterostructure from Fig. 1 for 
different values of annealing time (curves 2-4) and 
idealized step-wise approximation (curve 1). 
Increasing of number of curve corresponds to 
increasing of annealing time 
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Fig. 5a. Dimensionless optimal annealing time of 
infused dopant as a function of _ several 
parameters. Curve 1 describes the dependence of 
the annealing time on the relation a/L and = y=0 
for equal to each other values of dopant diffusion 


2167 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


coefficient in all parts of heterostructure. Curve 2 
describes the dependence of the annealing time on 
value of parameter ¢ for a/L=1/2 and €= y= 0. 
Curve 3 describes the dependence of the annealing 
time on value of parameter for a/L=1/2 and e= y 
= 0. Curve 4 describes the dependence of the 
annealing time on value of parameter y for 
a/L=1/2 and e= €=0 


We show optimal values of annealing time as 
functions of parameters on Figs. 5. It is known, 
that standard step of manufactured ion-doped 
structures is annealing of radiation defects. In the 
ideal case after finishing the annealing dopant 
achieves interface between layers of 
heterostructure. If the dopant has no enough time 
to achieve the interface, it is practicably to anneal 
the dopant additionally. The Fig. 5b shows the 
described dependences of optimal values of 
additional annealing time for the same parameters 
as for Fig. 5a. Necessity to anneal radiation defects 
leads to smaller values of optimal annealing of 
implanted dopant in comparison with optimal 
annealing time of infused dopant. 
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Fig.5b. Dimensionless optimal annealing time of 
implanted dopant as a function of several 
parameters. Curve 1 describes the dependence of 
the annealing time on the relation a/L and = y= 0 
for equal to each other values of dopant diffusion 
coefficient in all parts of heterostructure. Curve 2 
describes the dependence of the annealing time on 
value of parameter ¢ for a/L=1/2 and €= y= 0. 
Curve 3 describes the dependence of the annealing 
time on value of parameter for a/L=1/2 and ¢= y 
= 0. Curve 4 describes the dependence of the 
annealing time on value of parameter y for 
a/L=1/2 and e= €=0 


CONCLUSION 


In this paper, we introduce an approach to 
increase integration rate of element of a 
conventional static comparator. The approach 
gives us possibility to decrease area of the 
elements with smaller increasing of the element’s 
thickness. 
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Equations for the functions na 1), —, 3) and Vi. (y, 1; —, 3), i>0,j=0, k>0 and conditions for them 
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